A one-dimensional lattice model with exact matrix-product ground state 
describing the Laughlin wave function 
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We introduce one-dimensional lattice models with exact matrix-product ground states describing 
the fractional quantum Hall (FQH) states in Laughlin series (given by filling factors v = 1/q) on 
torus geometry. Surprisingly, the exactly solvable Hamiltonian has the same mathematical structure 
as that of the pseudo potential for the Laughlin wave function, and naturally derives the general 
properties of the Laughlin wave function such as the Zi properties of the FQH states and the 
fermion-boson relation. The obtained exact ground states have high overlaps with the Laughlin 
states and well describe their properties. Using the matrix product method, density functions and 
correlation functions are calculated analytically. Especially, obtained entanglement spectra reflects 
gapless edge states as was discussed by Li and Haldane. 

PACS numbers: 71.10.Pm, 75.10.Kt, 73.43.Cd 



I. INTRODUCTION 

The fractional quantum Hall (FQH) effect has been one 
of the most intriguing research field in condensed matter 
physics during the past few decades, in both theory and 
experiment Recently the FQH effect attracts much at- 
tentions due to new developments including the observa- 
tion in graphene^ flat band Chern insulators^, rapidly 
rotating bosons^ applications to topological quantum 
computing;* and so on. 

A key property of FQH states is their topological order. 
One consequence thereof is that their physical properties 
are insensitive to smooth deformations of the manifold on 
which we choose to study them. Using toroidal boundary 
conditions, the two-dimensional (2D) continuum system 
in a magnetic field can be reduced to a one-dimensional 
(ID) lattice model. This fact has been exploited in a se- 
ries of recent studies of the interacting many-body prob- 
lem in the limit geometry of a thin torus, referred as the 
Tao-Thouless (TT) limit 

In this paper, based on the ID approach, we intro- 
duce minimal models with exact ground states that de- 
scribe the FQH states in the Laughlin series (filling fac- 
tors v = 1/q). The physical quantities such as density 
functions and correlation functions are calculated analyt- 
ically by the matrix product (MP) method. The simplest 
case of our model (q = 3) reproduces a model describing 
the v = 1/3 FQH state that has been discussed quite 
recently by the authors with Bergholtaii and JansenJ^ 
We will also show that the structure of the ID mod- 
els [Eq. (|Sa)) below] implies that of the pseudo potential 
for the Laughlin wave function [Eq. (|10p below] , so that 
properties of the Laughlin wave function such as the Zi 
classification (the FQH effect is limited only odd q for 
fcrmions and even q for bosons) and general relations 
between fermionic and bosonic systems are naturally de- 
rived from our model. 

Among physical quantities obtained by the MP 
method, entanglement spectra (ES) attract the atten- 
tion recently. ES are obtained as coefficients when the 



wave function is divided into two subsystems. According 
to Li and HaldaneJ^, ES describe the "edge states" of 
the subsystem. Based on this notion, we extract (chiral) 
Tomonaga-Luttinger liquid behavior of the edge state an- 
alytically. 

This paper is organized as follows: In Sec. [Til we dis- 
cuss how ID models with exact ground states are related 
to the v — 1/q Laughlin states. We also consider bosonic 
version of the model. In Sec. IIII1 we introduce MP repre- 
sentation of the exact ground state, and perform analyt- 
ical calculation of the physical quantities such as density 
functions correlation functions and ES. In Appendix, we 
present the detailed calculation of the matrix elements of 
the pseudo potential, and those of the physical quantities 
based on the MP method. 



II. MODELS WITH EXACT GROUND STATES 

A. FQH system on a torus 

As discussed in the preceeding works jt~— when we 
consider 2D fermion system in a magnetic field using 
torus boundary condition with lengths Li in xj-dircctions 
(I = 1, 2), and single particle wave function with the Lan- 
dau gauge A = BX2X1 for the lowest Landau level, the 
system is described by the following model on a ID dis- 
crete lattice on direction x%, 



* = t E 

i=l | m |<fc<JV s /2 



C i+m C i+k C i+m+k C i> 



(1) 



where cj (cj) creates (annihilates) a fermion at site i, 
hi = eje^, and N s — L\L2/2-k is the number of sites (the 
magnetic length = \Jh/eB is set to be unity). The 
matrix element Vfc m specifies the amplitude for a process 
where particles with separation k + m hop m steps to 
a separation fc — to. This model conserves the center of 
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mass coordinate defined by K = X^=i jfij (mod N s ), 
which is related to the momentum of Xi direction. 

B. i/ = l/3 

In the preceding paper , 11 ! 12 a ID model with an ex- 
act ground state describing the v = 1/3 FQH state are 
constructed in the following way: First, as an approxi- 
mation in the vicinity of the TT limit, we truncate the 
Hamiltonian ([I]) at v = 1/3 up to the four leading terms 
(k + \m\ < 3). In the condition 

Vl = V 10 V 3Ql V 2Q > (2) 

which is satisfied for the v = 1/3 Laughlin state for any 
L\ with L2 — > oo, this model can be rewritten as 

Hi, a = , £,[Q\Qi + P}Pi\ (3) 

i 

where 

Qi = aic i+ ic i+2 + a CiC i+ 3, P t = l3 c % Ci + 2 (4) 

with ao,ai,/3o £ R. This Hamiltonian apparently has 
positive expectation values (Hi/ 3) > 0. Next, we intro- 
duce the wave function 

1*1/3) = II( ai " «° c li4 2 W*) l*?/3> > ( 5 ) 

i 

with |*° /3 ) = 1 100100100 •• •). Then the relation 

Qi 1*1/3) = Pi l*i/3) = 0,Vi is satisfied, so that |*i/ 3 ) 
is shown to be an exact ground state of Hi /a- This state 
is three-fold degenerate due to the center of mass conser- 
vation. This argument can also be applied to a v = 1/2 
bosonic FQH state with simple replacements. 

C. v = l/q 

Now we discuss how an exactly solvable model is con- 
structed for general Laughlin series v = 1/q with odd q. 
For v — 1/q with q > 3, the Hamiltonian should include 
longer range interactions than H1/3. The dimension of 
the subspace increases as (q+l)/2 = [q/2] + l where [■ ••] 
is the integer part of ■ • • (e.g. {|1001) , |0110)} for q = 3 
and {|100001) , |010010) , |001100)} for q = 5). We find 
that the extended model should be given in the following 
form, 

[9/2] 

1=1 i 
[9/2] 

Of = a^Ci+M^-M. ( 6b ) 

[9/2] -1 

p i l) ^ E pPd+^i+a-v-i- (6c) 

/i=0 



This is also positive semidefinite (Hi/ q ) > 0. We can 
construct a similar Hamiltonian with zero energy eigen- 
states if we take onl y sing le-/ term in (pa|) into account 
as was done in Refs. Il4lll5l . However, if there are not the 
degrees of freedoms of I, we cannot fix the parameters 
0$ uniquely, as will be discussed below. 

The corresponding wave function which gives the zero 
energy eigenstate is 

[9/2] 

i%) = IIE^™wi^)' ( ? ) 

i j=0 

where the squeezing operators are defined by Uk,o = 1 

and Uk,m(i) = for m + 0. This opera- 

tor is identical with a part of the pair hopping terms in 
Eq. (TT]). The root wave function |*?/ g ) is given by the 
charge-density- wave state in the TT limit, 

l*i /9 ) = l---i£ v _oiO v _o---), (8) 

9 9 

where every g-th sites is occupied. The parameters in 
Eqs. (|5a|) and ([7]) are determined so that they satisfy 
orthogonal conditions, oft> ■ s = 0, V7 where a' 1 ' = 

(«0 , a l J'" > a [q/2]) and S = C 1 ' S l>' " ' S [9/2])- SinCe 

the number of the subspace (components of s) is larger 
than the number that I may take by one, these equations 
fix the parameters s uniquely. On the other hand, there 
is no constraint for the parameter /3^\ since the oper- 
ators give pair hopping process between pairs sep- 
arated even number of lattice sites, which do not exist 
in Eq. 0. Thus Qf |¥jv_) = \* 1/q ) = 0,Vz is sat- 
isfied which means that ([7]) is the exact ground state of 
Eq. (|6a[) with zero energy. Due to the center of mass con- 
servation, the ground state is 9-fold degenerate, even if 
parameters a, [3 and s have site dependence. Uniqueness 
of the ground state in each center of mass sector can be 
shown considering that other states have finite positive 
energies. 

D. Bosonic systems 

The similar argument can also be applied to bosonic 
FQH systems with even q. The bosonic version of (UJ is 
given by 

H = E E V^UJi+Ji+m+kK ( 9 ) 

i=l \m\<k<N 3 /2 

where b\ (b/) creates (annihilates) a boson at site i. In 
this case the dimension of the subspace is (q + 2)/2 for 
even q (e.g. {|10001) , |01010) , |00200)} for q = 4). Our 
exact argument for fermions can be applied straightfor- 
wardly to the bosonic cases only replacing the fermionic 

operators by the bosonic ones q — > b 1 - in Eqs. © and 
©■ 
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E. Relationship with Laughlin states 

We discuss that our exactly solvable model well de- 
scribes the v = 1/q Laughlin states. In order to discuss 
the Laughlin sates in the present ID framework, we in- 
troduce a pseudo potential ; 16 ' 17 which is known to have 
the v = 1/q Laughlin state as the exact ground state, 



V(x) 



9-1 

E 

A=0 



c^V 2 ^ 2 ^), 



(10) 



where x = (xi,X2), c\ = (— 1) A |c> | are constants to keep 
the energy positive, and b > stands for the range of the 
interactions. Note that number of the terms increases as 
q is increased. This reminds us that our exactly solvable 
model (jfla)) has a similar structure. For N s L\, the 
single-particle wave function for the lowest Landau level 
in Landau gauge A = BX2X1 is given by 

where 2irk/Li is the momentum along the x\ -direction. 
In this basis, we derive the general form of the matrix el- 
ements of the Hamiltonian ([!} as the second quantization 

of m 



9-1 

vl =Vk' (a) - 

"km / j "km ' 
A=0 
A [(A-^)/2] 

53 E c»Ak 2 ~ m 2 y 

u=0 p,=0 



V, 



/(A) 



(12) 



km 



(k 2 +m 2 ) 



where are given by c\, 6, and L\ (see Appendix 
|A"| . Note that the function is antisymmetric (symmetric) 
under the exchange of k 2 and m 2 as a consequence of the 
Fermi (Bose) statistics. 

On the other hand, the matrix elements of our exactly 
solvable model Vk m (k + \m\ < q) are identified in the 
following way. First, we divide the matrix elements as 

Vkm = J2\ ^fem an d relate A to the variable in Eq. (paj) 
as A = 21 — 1 for fermions A = 21 for bosons. The matrix 
elements for each A which stem from [• • • ] in Eq. (|6aj) 
satisfy the relation 



(A) 



yW yW 
v k-m.O v k+m 



(13) 



This relation is a counterpart of Eq. ^ for the v = 1/3 
case, and is also satisfied if we choose the matrix elements 
in the following form 



9-1 

V km = E V kL' 

A=l,odd 
A— 0, even 



(14) 



V 



(A) 



--C x (k 2 



m 2 )V 



'-(k 2 +m 2 ) 
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FIG. 1: (color online) Overlap of the ground-state wave func- 
tion of the fermionic Laughlin state and that of our exactly 
solvable model J7J for q — 3,5 and JV e = N 3 /q as functions of 
L\ in finite-size systems, calculated by exact diagonalization. 
The data are plotted up to the isotropic points of the torus 
(Z/i = L%). The approximation becomes better for larger-g 
cases, since the number of long-range terms are added as q is 
increased. 



where C\ is a constant and A is restricted to be odd 
(even) for fermions (bosons) due to the statistics, as al- 
ready mentioned. Although Eqs. (|12l) and (p~4|) coincide 
at q = 3 (q = 2) for fermions (bosons), they are differ- 
ent in general. However, we may interpret the truncated 
Hamiltonian of ([1]) with (|14[) as an approximation of the 
full Hamiltonian with (|12l) for L\ — » limit (see Ap- 
pendix |Aj . 

In order to show that our model well describes the 
Laughlin state, we numerically calculate overlap of the 
ground-state wave function of the Laughlin state and 
that of our exactly solvable model Q in finite-size sys- 
tems using exact diagonalization. As shown in Fig.[T] the 
overlap stays close to 1 up to L\ ~ 6 for all cases, and 
the high overlap regions are extended for larger q. Thus 
we may conclude that our wave function well describes 
the Laughlin states at v = 1/q with q > 3 (q > 2) for 
fermions (bosons). Note that the approximation becomes 
better for larger-g cases, since the number of long-range 
terms are added as q is increased. This is in contrast with 
the AKLT model as an approximation of the integer-spin 
Heisenbcrg model which includes more extra terms as the 
length of the spin is increased. 



F. Z-2 properties 

The present exact result seems to be applicable to fill- 
ings v = 1/q with even (odd) denominator q for fermions 
(bosons) as well. However, the structure of our model 
naturally excludes such situation, in accordance with the 
symmetry of the Laughlin wave function! 2 - For fermions 
(bosons) with even (odd) q, the expressions of Qf \ P^ l \ 



and Eq. (|14[) remain unchanged, while the summations 
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in Eqs. (|6a[) and (0 are replaced by Ya=i ^ an d 
Y^jlo 2 \ respectively. This is because the number of 
subspace for fermions (bosons) with even (odd) q system 
becomes [(q + l)/2] (e.g. {| 10001) , |01010)} for q = 4 
fermions). Thus the number that I may take becomes the 
same as the number of the subspace. Therefore, the pa- 
rameter s cannot be determined uniquely for these cases 
by the equation = 0,Vi, which means that models 

with the exact-ground states (0 cannot be constructed. 



G. Compressibility 

We can show the vanishing of the compressibility in 
this system as is expected for FQH states. As discussed 
in Ref. [Til, zero energy eigenstate can be created insert- 
ing everywhere in \^i/ q )- This process makes a par- 
ticle pair separated q + 1 sites while the range of in- 
teractions is q sites, so that it is equivalent to making 
an open boundary. Therefore the ground-state energy 
of such cases gives E q (N s + 1) = E q (N s ) = 0. On the 
other hand, if we shrink the system size by removing 0, 
a defect corresponding to the fractional charge e* = e/q 
appears. Since this defect localizes, the excitation energy 
E q (N s — 1) has very small size dependence as was explic- 
itly shown in the v = 1/3 case^i Therefore, the inverse 
compressibility 



1 N 

- = lim -j T [E(N s -l)+E(N s + l)-2E s (N s )} (15) 



clearly diverges. 



III. MATRIX PRODUCT METHOD 

From the exact solution ([7]) it is possible to calculate 
generic ground-state expectation values. For this pur- 
pose, it is convenient to introduce a matrix product (MP) 



represcntatio: 



,18,19 



of the ground-state wave function for 



v = 1/q FQH states with q = 2m + 1 (q = 2m) for 
fermions (bosons). For a periodic system with N s = qN 
sites, ([7]) can be written as 




FIG. 2: (color online) (a) Density functions {fik) and (b) cor- 
relation lengths £ of the density-density correlation functions 
of v — 1/q FQH states with q — 3, 5, 7 in L2 — s- 00 limit 
as functions of the circumference of the torus L\, calculated 
analytically by the MP method. 



Here we have introduced the spin-m representation for 
g-sites unit cell for fermions (bosons), 



|10 • • ■ ■ 
|01 - ■ • 
|00 ■ - ■ 1 • 
loo - • • 



■00) (|00 ■ 

■00) ([10 • 

■00) (|00 • 

■01) (|00 • 



•00---00)) 
■00- ■■00)) 
•10---00)) 
•00---02)) 



(18) 



The matrix (fTT)) has the same size as that of the valence- 
bond-solid states of an S — m quantum spin chain^S 
however, there isaraxm zero entry part and the transfer 
matrix G(ri,r 2 ;t 1 ,t 2 ) ~ 
eigenvalues A± 



9j{ri,ti)gj{r2,t 2 ) has only two 



(l±-\/4s 2 — 3)/2. Due to this property, 
the matrix G can be reduced to 4 x 4 form for all q, so 
that the calculation of expectation values is extremely 
simplified. 



Correlation functions 



= tr [5i52 •••9n], 



(16) 



where gi is identified as the following (m + 1) x (m + 1) 
matrix (see Appendix |B|) . 



*i|-l>< 
«a|-2>< 



-1), 1+2) 





-m) 





(17) 



Using the MP formalism we obtain the density func- 
tion in the infinite systems N — > 00, which has g-site 
periodicity as (see Appendix IC 2j) . 



2(s 2 - 1) 
1 

V4s 2 - 3' 



1 



1 



V4l 



(19) 
(20) 



Figure[5](a) shows the L\ dependence of the density func- 
tions for q = 3, 5, 7. Since our wave function is not trans- 
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<?= 5 



torus • cylinder 



?= 3 



-1 1 



-1 1 




FIG. 3: (color online) The entanglement spectra of the 
fcrmion systems (q — 3, 5, 7, N — > oo) on torus (cylinder) 
geometries at Li = 9 as functions of the momentum for 
ii-direction AK. They form the conformal tower structure 
of (chiral) Tomonaga-Luttinger liquids describing the edge 
states. 



function ((7]) is written as the Schmidt decomposition 
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-* /2 N&>®rV&> 



(22) 



[<5f 5 ]- 1/2 [ntift( + L)ki a are or- 



where IV^^) 

thogonal states in the subsystem -A(-B), and A/^ J2 
are their norms. The ES are given by £j = 
— \og{j\fijj\ff 1 /j\f). As was pointed out by Li and 
Haldane^ the ES characterize the "edge states" of the 
subsystem. Although the number of the ES is limited as 
m + 1, the number of entanglement levels for bulk FQH 
states grows with the size of the systemJ^ However, when 
the ES arc plotted as functions of the center of mass AK 
(which corresponds to momentum of X\ direction) of the 
subsystem as shown in Fig. [31 they give low-lying liner 
spectra of the gapless edge state described as a chiral 
Tomonaga-Luttinger liquid.— 

Next, we evaluate ES in the torus geometry where the 
Schmidt decomposition is given by 



Mi 



(23) 



lationally invariant, it does not describe a genuine liquid 
state, however, these results indicate that the systems ap- 
proach to homogeneous liquid states as L\ is increased. 
This is consistent with the result obtained by Rezayi and 
Haldane which showed that the charge-ordered state in 
the TT limit is adiabatically connected to the FQH liquid 
state through the smooth deformation of the torus4 Ac- 
tually, our exact result describes many important aspects 
of the FQH states, as we have already seen. Similarly, the 
density-density correlation functions are calculated which 
have exponentially decaying forms as (see Appendix IC 2[) 



1 - V4s 2 - 3 



(21) 

where £ is the correlation length. As shown in Fig. EJb), 
the correlation length becomes longer in large L\ regions 
as q is increased. Since the correlation length is pro- 
portional to the inverse of the energy gap, this result is 
consistent with the fact that the v = 1/q FQH state be- 
comes unstable for larger q. Similar behavior has also 
been confirmed in v = p/(2p + 1) Jain series^ 



B. Entanglement spectra 

The MP formalism enables us to derive entanglement 
spectra (ES) of the FQH system {£j} analytically (see 
Appendix [DJ ■ First we consider the cylinder geometry 
(open boundary conditions), and divide the system into 
two blocks A and B with length L = N/2. Then the wave 



The ES &j = - \og{N^N® /N) are characterized by two 
separated "edge states" of the subsystems. The number 
of the ES is (m + l) 2 , and m(m + l)/2 of them are two- 
fold degenerate. As shown in Fig. [3l the ES show a "di- 
amond shape" conformal tower structure of a non-chiral 
Tomonaga-Luttinger liquid which has both left and right 
movers^ This result can be explained by the relation 



obtained by the MP method 

((AA'i - AK 2 ) = £(Aifx) + £(AJr a ) - £o 



(24) 



for TV — > oo, where £o is the lowest ES. This means that 
the ES are given by the spectra of an effective noninter- 
acting ID "entanglement Hamiltonian" with right and 
left modes with momenta AKi, ~AK 2 . We can also 




FIG. 4: (color online) The entanglement spectra and entan- 
glement entropy of the fermion systems (q = 5, N — > oo) on 
torus geometry as functions of the circumference of the torus. 
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evaluate the difference between ES on a torus and those 
of cylinder as log[A + /(A + — A_)]. 

In Fig. HI L\ dependence of the ES and entanglement 
entropy (EE) for q = 5 with torus boundary conditions 
are shown. The diamond structure of ES is broken for the 
large L\ regions and the EE S approach to log(m + l) 2 . 

IV. CONCLUSION 

We have introduced a ID minimal reference model with 
exact ground states which describes v = 1/q FQH states 
with odd (even) q for fcrmions (bosons). This model is 
constructed as the second quantization of the pseudo po- 
tential for the Laughlin wave function, around the TT 
limit of the torus boundary system. The model has the 



same structures of the pseudo potential, and naturally 
derives general properties of the Laughlin wave function. 
The ground states are described by the MP wave func- 
tions that enable us to calculate physical quantities an- 
alytically. Especially, we have analytically obtained the 
conformal tower structure of the ES which reflects chiral 
Tomonaga-Luttinger liquid behavior of the edge state. 
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Appendix A: Matrix elements of the pseudo potential Eqs. (]12[) and (|14[) 



We consider a model of N interacting electrons on a torus with lengths Li in ^/-directions (I = 1,2). When the torus 
is pierced by N s magnetic flux quanta, boundary conditions requires LiL 2 = 2ttN s . In Landau gauge, A = Bx 2 x\, a 
complete basis of N s degenerate single-particle states in the lowest Landau level, labeled by k = 0, . . . , N s — 1, can be 
chosen as 



Eexp 



m£2 



/ 27F , r 

1 [ + L 2 m \ X\ 



exp 



27r , r 
7; ( x 2 + —k + L 2 m 



(Al) 



where unit of xi and Li is the magnetic length ^b- In L-2 — > 00 limit, it becomes Eq. In this basis, any 

translation-invariant 2D Hamiltonian with two-body interactions assumes the following ID lattice model, 



fcifc 2 fc 3 fc 4 |m|</c<7V s /2 



Vkm — Vkm C i+m C i+k C i+k+m C i 

i 

= £ V i+m , i+k , i+m+k>i T <Whww» + V i+k , i+m ,i,i +m+k T Vi +k , i+m , i+m+k ,i, 



(A2) 



where the upper (lower) signs correspond to fermions (bosons) and the operators c k (c\.) annihilates (creates) a particle 
at site (orbital) k. The matrix elements V qi q 2q3qi can be calculated according to the standard second quantization 

1 



y fcl fe 2 fe 3 fe 4 = Vk^k s kJ h J h7 c. k c u = ^JJ ?Afc 1 (a;i)V4 2 ( a; 2)V A (a;i - x 2 )^ k3 (x 2 )tp ki (x^x^x?, (A3) 

where ^ k {x) = i[> k (x)c k . One can find that for fermionic (bosonic) systems the matrix elements are real and antisym- 
metric (symmetric) functions for k and m reflecting Fermi (Bose) statistics. Using the Fourier transform V{p) as the 
interaction potential V{x\ — x 2 ), the matrix elements Vk 1 k 2 k 3 k i are given by 



*^k\k 2 k^ k& 



J ki+k 2 ,k 3 +k 



2LiL 2 



± Y. s 'ki-k i , P1 L 1 /2^(py 



p 2 /2-i(k 1 -k 3 )p 2 L 2 /N a 



(A4) 



where p = t (pi,p 2 ), pi = 2nmi/ L±, p 2 — 2i:m 2 / L 2 , mi,m 2 £ Z, and 8' is the Kronecker delta function with a 
period N s . The pseudo potential proposed by Trugman and Kivelson is known that it has the v = 1/q Laughlin wave 
function as the exact ground stated 



q-l q-1 

V(x) = V (A) = Y cxb 2X V 2X 5 2 {x), x = x x - x 2 , 

A=0 A=0 



(A5) 
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where b is the range of the interaction. For N s ^> L±, the matrix elements Vkm are calculated in the L2 — > 00 limit 
using the binomial theorem and Gauss integrals 5 n = x 2n e~ x I 2 as 



V 



(A) _ {-b^Yc\ _ (K 2 +M : 



)/2 



47rii 
x dx 



(M + x ) (e"^ 

A 



2 -I- t- 2n I A ( p-(^-^) /2 j. p -(x+iK) 2 /2^ T ^2 + ^A f p -(x-iM) 2 /2 _,_ p -(.T+iM) 2 /2 



)/2 



47rLi 



x dx 



J2 \C U (K 2 -M 2 y 



v=0 



(-iy(K 2 + X 2 ) X -» ( e -( X " iK ) 2 / 2 + e"( a: + iif ) 2 /2^ zp ( M 2 + ^jA-v / e -(x-iM) 2 /2 + e -(x+iM) 2 /: 



: ( Z ^£A e -(, 
47rii 



(~ b ) C A „_(g2 +M 2 



)/2 



2V27rii 



[(A-i/)/2] 

! Y, aC,(K 2 - Af 2 )" (-l)^ 2 ^ 1 a _ v C 2m 5 A - v -^ [(-l) 1 ^ 2 " T M 2 ^] 

!/=0 fi=0 

A [(A-f)/2] 

y y (_i)^ 2 2m+i ^ ; y tt L ± < K * _ M 2 r rr-i)^ T m 2 1 



jy = (U=0 



where if = 2 pp and M = ^j^ 1 - One can find that in fermion (boson) systems the powers of (K 2 — M 2 ) are always 
odd (even) for each terms, and the leading term of K . M appears as (K 2 — M 2 ) A when A is odd (even), because 
[(— X) V K 2 ^ =F M 2/i ] vanishes for (y,fi) — (A,0) when A is even (odd). Therefore Eq. (|12l) can be approximated as 
Eq. (TT4")) in the small L\ regions. 



Appendix B: Matrix Product Representation 
(Eq. (9)) 



Let us consider the case oi v = 1/3 as the sim- 
plest example. The ground state manifold of this 
model is three-fold degenerate, spanned by charge or- 
dered states with one electron per a three-site unit cell: 
| - • • 010010010 •••). The V2i-term induces fluctuations 
upon these ground states through the process, 



|010 010) O |001 100) . 



(Bl) 



The truncated model can be mapped to an S — 1 quan- 
tum spin chain by identifying the states of the unit cell as 
|010) -> |o), |001) -> |+), and |100) ->■ |-). We consider 
the matrix product representation of the ground state 
wave function in the S = 1 spin basis. Then the matrix 
can be written in a following 3x3 form 



f—Hi f-o\0) t f-+\+)i 
fo-\-)i /oo|o) 4 /o+|+) 4 

/ + _|-). z+oio), / ++ |+>. 



(B2) 



where i is the site index for the unit cell and /io are con- 
stants. The index I and O denote the input and output 
states that are the configurations of the (i — l)-th and 
the i-th spins. We suppose that the ground-state wave 



function in periodic boundary conditions is given by 

N 

l*i/3> p =I[( 1 + 10100100100...) , (B3) 

j=i 

N 

= JJ(l + «iZ3i)|ooo...} 

i=l 

=t r [fflff2 ' --9N], 

where si = ao/ai, and 



(B4) 
(B5) 



^fc — C I+l C l+2 C fc+3 C fe' 



SfSi Si+iS i+1 . (B6) 



One can find that the states of the nearest two spins can 
be written as 

ifii-! ®tpi = |oo) + si |+-) + |o+) + |-o) + |-+) . (B7) 
Therefore, we can set 

: 



/ — Jo- — f+o — /+ + 

foo /o+ /— o / — h 
/+- = Si. 



1 



Thus, the matrix Eq. (|B2j) then becomes 

|o>, |+), 
lo), |+) 
»i|->, 



(B8) 
(B9) 
(B10) 



(Bll) 



If we change the matrix (|B2|) in different basis (s x |— ) . + 
> l+)t)j ^ ne above matrix can be reduced to 2 x 2 form 
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without losing generality, 



fj> 



si 







(B12) 



Although Eqs. (|B11|) and (|B12|) give the same result, we 
use the latter representation for simplicity. 

In the diagonal component of the matrix product, the 
upper one is both for open and periodic boundary con- 
ditions, while the lower one is only for periodic bound- 
ary systems. Therefore, for open boundary systems, the 
ground-state wave function is given as 



l*i/3> = tr'[.9i.g 2 • • -g N ], 



(B13) 



where tr' means partial trace for the upper component. 

For general v = 1/q with q = 2m + 1, we consider in a 
similar way as the above using the spin-m representation 
(|18|) . and obtain a q x q matrix as 









-to) 







|o) |+1) 

|o) |+1) 

-1) •■■ 

... o ••• 



\-m) 








(B14) 



The matrix with the reduced dimension to + 1 can be 
obtained as Eq. (JT7J) in the text. 



Appendix C: Correlation functions 

1. Transfer matrix 

We consider the matrix product representation of the 
v = l/(2m + 1) ground state wave function in (2to+ 1)7V 
sites periodic system. The normalized ground state wave 
function can be written as 

|*m) = W~ 1/2 Tr[gig 2 ■ ■ ■ 9n], (CI) 
jV=Tt[G n ], G = g t ®g t , (C2) 

where gi is identified as Eq. ([T7]). The (m + l) 2 x(m + l) 2 
transfer matrix G(n, r 2 ; t\, t 2 ) = g\{r\,ti)gj(r 2 ,t 2 ) has 
only two eigenvalues, 



We define (to + 2) x (to + 2) local matrix 1 to reduce the 
dimension of G to (m + 1) x (m + 1): 



G 



1 1 



(sfl) ••• 
(41) •■• 



1 








••• 



(C5) 



Thus the transfer matrix is diagonalizcd by finding a uni- 
tary matrix U as 



G = U^GU 



(A_l) ■•• 
(A+l) 







•• 











(C6) 



Thus we can calculate physical quantities just in the same 
way as in the case of v ~ 1/3.— 



2. Density and Charge correlation 

In the spin-m representation for the v = 1/(2to + 1) 
state (|18p. the density operator and its transfer matrix 
can be written as 



'■qi 



3 = S Sf 



St J' 



G[n](j) 



,S Srj3i 



(C7) 



Using the unitary transformation U defined in Eq. (|C6[) . 
we obtain the density function with g-sites periodicity as, 

V(j;N) EE {flqi+j) = (hqi-j) 

Tr [U-iG[n](\j\)UG»- 1 ] 



lim V(j;N) 

N— J-oo 



2(s a -l) 
1 



Tr[G» 
1 



Troy CjV°) 
(j = o) 



(C8) 



The charge-charge correlation function is calculated as, 
C(i,j,i',j';N) ee (h qi+j h qt , +r ) 



Tr 



Tr [Go 



ATI 



(C9) 



A = 0, 



A+ = 



(C3) 



that are obtained from the following secular equation, 

m 

= dot |G - A| = (1 - A)(-A) m - (-A)™" 1 s r (C4) 



lim C(i,j,i',f;N)~ 

N,\i—i'\-¥oo \A+ 



'h^l M Ee -#-i'i/s 

1 + \/4s 2 - 3 / 



(C10) 



where £ is the correlation length for the infinite system. 
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Appendix D: Entanglement Spectra and Entropy 

In this section, we compute the entanglement spectrum 
(ES) and the entanglement entropy (EE) both in cylinder 
and torus geometries. We obtain the conformal tower 
structure of the ES in Fig. [3] analytically. 



1. Cylinder 



where J\f^ B ^ is the normalization factors that is given 
by 

Mf ee [G'G^ 1 ] j(m+2) _ (ro+1)jl , (D10) 
X?=[G N - L - 2 G'\ j[m+2) _ {m+iy (Dll) 
The reduced density matrix is given by 
M A M B 



For a open boundary system, the normalized ground 
state wave function can be written as 

|tt m ) =Af- 1/2 Tr[g[ 92 ---g'^ J\f = Tr [G'G N ~ 2 G"] , 

(Dl) 

where the terminators g[,g'{ are identified as the follow- 
ing matrices 



|0> i 












|0> 4 

Si 1-1) 



-m) 




s m \-m) i 



and the transfer matrices G', G" are given by 



G'=g^g' = 



= 9i 



1 1 




■•• 

1 ••• 

(sfl) ••• 

(4,1) o ••• o 



(D2) 



(D3) 



(D4) 



(D5) 



We cut this system into blocks A and B with length L 
and N — L, respectively. The wave function of the whole 
system is divided as 



The ES if and EE S A = -Tr A [p A log p A \ can be written 



as 



S J 



( AffjVf \ 



(D13) 
(D14) 



The normalization factor for the wave function is calcu- 
lated by using Eq. (|D1[) as 



N = 



A_ - A i ' 



(D15) 



For L = N/2 => G L = G N ~ L = G N / 2 , the normalization 



factors Afj , Afj are written as 



X L + 1_ X L + 1 

A_-A + 
A--A+ 

A L + 1_ A L + 1 

A_-A + 



3 1 X -~ X + -2 

\--\+ "j 



U = 1) 
(J > 1) 

(J = 1) 

V > !) 



(D16) 



(D17) 



We now can obtain the ES and EE by substituting the 
j\ff,j\ff into Eq. (|DT3jl . For L = N/2 -> oo, 



A_-Aj 



(X.-X + )X + °3 



(J = 1) 

-i (i>i) 



(D18) 



When the ES are plotted as functions of the shift of the 
center of mass AK of the subsystem, it shows approxi- 
mately a liner relation which corresponds to a conformal 
tower of the chiral Tomonaga-Luttinger liquid. 



\* m }=Ar- 1 /*Tr\ g ' l92 .--g't f ] 



m+l 



(D6) 



: N V2 H {[3i52 ' [9L+19L+2 ■ ••'/",.,} 



We can then define the block states as 



\^f)^{Mf)-^[g' 1 g,---g L \ 1 ., 

\^!)^{NZ)-^[g L+l g L+2 ...g^} 



(D7) 



(D8) 
(D9) 



Torus 



We cut this system into blocks A and B with lengths 
L and N — L, respectively. The wave function (|C2j) of 
the whole system can be divided as 



*m) =M- 1/2 Ti[ gi g 2 ---g N ] 



rn+1 



AT 1 / 2 £ {[g l92 ---g L ] 



9\9i ■ ■ ■ 9L\ jlj2 L3i+i5i+2 • • • gN\j 



31 ,32=1 



(D19) 

2jl } ' 

(D20) 
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We can then define the block states as 

\^ = WLr 1/2 &92---9L] jlj2 , (D21) 
fe) = (Mf lh )-V 2 [ 9L+l9L+2 --- 9N ]. 2ji , (D22) 

where jVj^f ^ s the normalization factors for A(B) block, 

■^jlj-2 = [ G ]j 1 ( m+ 2)_(m+l)j 2 (m+2)-(m+l) ' (D23) 
= ] j 2 ( m+ 2)-( m +l)j 1 ( m +2)-(m+l) ' (^24) 

One can now rewrite the ground state wave function as 
Schmidt decomposition, 



Jlj2 V 



J2* JU2 1 I / A \ ^ |„/,B 



TV 



M& a >®M& a >. (D25) 



The reduced density matrix is given by 

The ES f ^ j and EE S A can be written as 
£ A ■ = - log nn nn 



(D26) 



(D27) 



AT A - AT B ■ ( \f A Af B ■ \ 
S A = J2 J1J f/ lj2 log "m"™ , (D28) 



3132 



Let us consider the normalization factors j\f and 
r A(B) r , L 

'3132 



Aff^fl. G L in Eq. (|D23|) can be calculated as Therefore, 



for L = N/2 =>■ G L = G W ~ L = G w / 2 , the normalization 
factors M A j2 , 7V^ J2 are given by 



j\f A . = 



3l32 



jV B ■ 

J 3132 



A_-A+ j'i-l 

A_-A + 
A^-A^ 1 



A_-A 4 



(ii = i, h = i) 
(ii > i, .?2 = i) 

0i = 1, .72 > 1) 

(il > 1, 32 > 1) 



: j\f A - . 
J 3231 



(D29) 



(D30) 



Now we can use Eq. (|D27[) to calculate the ES and the 
EE. For L = N/2 -> oo, we get 



M A - 

■ 3132 



77 



A 








A 








A 


-A+ 




a; 1 


A- 





(ii = i, .72 = i) 

(ii > i, 32 = i) 

(ii = 1, .72 > 1) 

(ii > i, 32 > i) 



3132 J 3231 



(D31) 



(D32) 



When the ES are plotted as functions of the shift of the 
center of mass AK of the subsystem, for N — > oo, we 
can show the relation (|24[) which means that the ES of 
the torus are given by the combination of the right and 
the left movers of the chiral Tomonaga-Luttinger liquids. 
Therefore, the ES shows the conformal tower structure 
of the non-chiral Tomonaga-Luttinger liquid. We can 
estimate the difference of the ES on a cylinder and that 
on a torus as log[A+/(A+ — A-)]. 
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